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ABSTRACT 


This  study  is  concerned  with  the  problem  of  calculating 
potential  flow  about  a  nonlifting  body  in  an  unbounded  fluid. 
Several  simple  explicit  approximations  for  the  velocity  potential 
are  obtained  and  investigated  numerically.  Results  of  calculations 
are  presented  for  the  simple  cases  of  potential  flows  due  to 
translations  of  ellipsoids  and  ogives. 
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1.  Introduction 

A  classical  and  fundamental  problem  in  hydrodynamics  is  that  of  calcu¬ 
lating  potential  flow  due  to  motion  of  an  arbitrary  body  through  an  unbounded 
ideal  (incompressible  and  inviscid)  fluid.  This  potential-flow  problem  in¬ 
deed  is  directly  relevant  to  naval  hydrodynamics  (e.g.  calculation  of  flow 
about  a  deeply-submerged  body  and  about  a  ship  in  slow  motion)  and  to  subsonic 
aerodynamics.  Mathematically,  the  problem  consists  in  determining  the  velocity 
potential  ($)  given  by  the  solution  of  the  Laplace  equation  7  $=*0  in  the  flow 
domain  (d)  exterior  to  the  body  surface  (b) ,  subject  to  the  condition  that 
<|>  vanishes  "at  infinity"  (far  away  from  the  body)  and  to  the  Neumann  condition 
that  the  normal  derivative  <)>n  of  the  potential  is  prescribed  on  the  body  sur¬ 
face. 

A  great  variety  of  methods  have  been  developed  for  solving  the  above- 
defined  exterior  Neumann  boundary  value  problem.  These  methods  can  be  divided 
into  two  basic  classes  of  methods,  namely  integral-equation  methods  in  which 
an  integral  equation  is  formulated  and  solved  numerically,  and  direct  numerical 
methods,  mainly  of  the  finite-element  type;  hybrid  methods,  in  which  finite 
elements  and  an  integral  representation  are  coupled,  however  have  also  been 
used,  e.g.  by  Jam!  and  Lenoir  (1977).  Integral-equation  methods  are  widely 
favored,  as  is  apparent  from  the  long  list  of  references  in  the  review  of  Hess 
(1975)  and  Kcfmer  and  Hirschel  (1977),  and  indeed  are  natural  and  a-priori 
advantageous  methods  in  the  sense  that  they  Explicitly  take  advantage  of  the 
Laplace  equation  for  reducing  the  original  three-dimensional  problem  to  a  two- 
dimensional  computational  problem  in  which  the  velocity  potential  only  needs 
to  be  calculated  on  the  body  surface. 

One  of  the  earliest  and  most-widely  used  integral -equation  methods  is 
that  of  Hess  and  Smith  (1966) ,  in  which  an  auxiliary  distribution  of  sources 
on  the  body  surface  is  employed,  in  the  classical  manner  described  by  Kellogg 
(1929).  Modifications  of  this  method  have  been  devised  by  Landweber  and 
Macagno  (1969),  whose  method  mainly  differs  from  that  of  Hess  and  Smith  in 
the  treatment  of  the  singularity  of  the  kernel  of  the  integral  equation  and 
in  the  procedure  for  obtaining  numerical  solutions,  and  by  Webster  (1975), 
whose  method  is  based  on  a  Fredholm  integral  equation  of  the  first  kind  ob¬ 
tained  by  assuming  a  distribution  of  sources  on  a  surface  slightly  inside 
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the  body  surface.  Methods  based  on  the  use  of  an  auxiliary  distribution  of 
normal  dipoles  (instead  of  sources)  on  the  body  surface  have  also  been  used, 
e.g.  by  Chang  and  Pien  (1975).  Finally,  methods  based  on  an  integral  equation 
for  the  velocity  potential  itself  (rather  than  for  assumed  auxiliary  distri¬ 
butions  of  sources  or  dipoles  as  in  the  previous  methods)  have  recently  been 
used  by  Chow,  Hou,  and  Landweber  (1976),  and  by  Noblesse  for  the  analogous 
but  more  complex  problems  of  ship  wave  resistance  (1978)  and  of  flow  about 
a  body  in  regular  water  waves  (1980) . 

In  this  study,  a  new  basic  integral  identity  is  presented.  An  inter¬ 
esting  feature  of  this  identity  is  that  it  is  valid  in  the  whole  space  (that 
is,  both  in  the  flow  domain  and  inside  the  body,  as  well  as  exactly  on  the 
body  surface) ,  and  is  in  fact  equivalent  to  a  set  of  three  classical  integral 
identities  that  are  exclusively  valid  for  points  outside  the  body,  exactly 
on  the  body,  and  inside  the  body.  In  the  particular  case  when  the  Laplace 
equation  is  satisfied,  this  "unified  integral  identity"  becomes  identical  to 
the  integral  equation  of  Chow,  Hou,  and  Landweber  (19^6).  This  integral  equation 
is  first  investigated  by  considering  the  simple  cases  of  potential  flows  due  to 
translations  of  elliptical  cylinders  and  ellipsoids.  The  investigation  then 
is  continued  for  other,  complementary,  simple  cases,  namely  that  of  potential 
flows  due  to  translations  of  ogives.  In  particular,  several  simple  explicit 
approximations  for  the  velocity  potential,  which  are  suggested  by  the  above- 
mentioned  study  of  flows  about  ellipses  and  ellipsoids  and  indeed  are  exact  for 
these  particular  flows,  are  evaluated  numerically  and  compared  to  the  exact 
potential  (given  analytically  by  conformal  mapping).  Finally,  the  basic 
integral  equation  of  Landweber  is  expressed  in  a  «  jdified  form,  also  of  the  non- 
homogeneous  Fredholm  type.  The  given  term  in  this  modified  nonhomogeneous  Fredholm 
integral  equation  provides  another  explicit  approximation  for  the  potential.  This 
explicit  approximation  is  also  exact  in  the  particular  cases  of  translations  of 
ellipsoids,  and  it  is  compared  to  the  exact  potential  for  translations  of  ogives. 
The  calculations  for  translations  of  ogives  show  that  fairly-good  approximations 
to  the  velocity  potential  can  be  obtained  by  the  above-mentioned  simple  explicit 
approximations,  notably  the  last  one.  For  practical  purposes,  the  main  results 
of  the  present  study  indeed  reside  in  these  simple  explicit  approximations. 
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2.  Basic  integral  identities 

Througtout  this  study,  variables  are  supposed  to  be  nondimensional. 

It  is  well  known  that  fundamental  integral  identities  can  be  obtained  for  the 
velocity  potential  by  applying  a  classical  Green  identity  to  the  potential 
function  4> (x)  and  the  free-space  Green  function  G(x,£),  which  is  given  by 

G(x,|)  =  -l/4*[(x-02+(y-n)2+(z-02]1/2  •  (2.1) 

This  classical  Green  identity  is 


[  (<t.V2G-GV24>)dv  =  \  (4>G  -G<f>  )da,  (2.2) 

U  ^b  n  n 

222  1/2 

where  the  fact  that  <j>=0(l/r)  and  G=0(l/r)  as  ri(x  +y  +z  )  was 

3 

used  to  discard  the  integral  of  <pG  -G$  =0(l/r  )  over  a  large  surrounding 

n  £  n 

sphere  of  radius  r  (surface  area  ^  r  ).  In  formula  (2.2),  and  indeed  here- 
after  in  this  study,  the  following  notation  is  used:  (x) ,  G=G(x,£), 

v=  (3x»3y»8z^  ’  Gn-3G/3nHVG‘a  and  4>n=34>/9nHV(fcn,  where  n=n(x)  is  the  unit 
normal  vector,  at  point  x,  to  the  body  surface  b  pointing  inside  the  body 
(outside  the  flow  domain  d) ,  as  is  shown  in  figure  1;  furthermore,  dvidv(x) 
and  daida(x)  represent  the  differential  elements  of  volume  and  area  at 
point  x  of  d  and  b,  respectively. 

The  Green  function  G(x,C)  satisfies  the  well-known  equation 


s2g  =  6(x-C)6(y-n)6(z-0 


(2.3) 


so  that  we  have  4>  (x)V2G=<}>  (t)V2G+[4>  (x)H>  (f )  ]V2Gi$  (t )V2G,  provided  $ 
is  continuous  in  d+b,  which  is  assumed  here.  Equation  (2.3)  then  becomes 


C4>* 


GV2<|>dv 

d 


(G$  — <J> G  )da 
n  n 


(2.4) 


where  the  notation  4>A=$(C)  was  used  for  shortness  (and  will  indeed  be 
used  hereafter  in  this  study),  and  C  is  defined  as 
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(2.4a) 


C 


V2Gdv  . 
d 


It  follows  from  equation  (2.3)  that  we  have  C=1  if  the  point  £  is  in  the  ex¬ 
terior  domain  d  strictly  outside  the  body  surface  b,  whereas  we  have  ChO  if 
£  is  in  the  interior  domain  d.^  strictly  inside  b. 

We  thus  obtain  the  classical  integral  identities . 


GV  $dv  - 


(G4)  — G  )da 
n  n 


(2.5a) 


for  £  in  d-b,  and 


GV  <f>dv  - 


(G<{>  — <£G  )da 
n  n 


(2.5b) 


for  £  in  d^-b.  It  can  also  be  seen  from  equation  (2.3)  that  we  have  01/2 
if  the  point  X  is  exactly  on  the  body  surface  b,  at  least  for  points  X 
where  the  surface  b  is  smooth;  more  generally,  the  value  of  4ttC  at  a  point 
£  of  b  is  equal  to  the  solid  angle  under  which  the  exterior  domain  d  is 
viewed  from  £.  We  then  have 

(G<,ln_((,Gn)da  (2.5c) 

for  £  exactly  on  (smooth)  b.  Identities  (2.5a)  and  (2.5c)  are  usually  obtained 
by  applying  the  Green  identity  (2.2)  to  the  domain  d-e,  where  e  represents 
the  domain  inside  a  small  sphere  centered  at  the  point  X,  or  the  inter¬ 
section  domain  of  this  sphere  with  d  if  X  is  exactly  on  the  boundary  sur- 

2 

face  b;  we  then  have  V  G50  in  d-e,  while  on  the  right  side  of  equation  (2) 
there  is  an  additional  integral  over  the  surface  of  the  small  sphere  (or 
portion  thereof)  whose  value  can  be  shown  to  be  equal  to  -$A(or  ~4>*/2  if 
X  is  on  b)  in  the  limit  as  the  radius  of  the  sphere  vanishes.  This  tradi¬ 
tional  derivation  and  the  derivation  used  in  this  study,  based  on  equation 
(2.3),  are  of  course  equivalent. 


2**  “  1, 


GV  4>dv  - 
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As  is  explicitly  indicated  in  equations  (2.5a,b,c),  the  value  of  the 

constant  C  on  the  left  side  of  equation  (2.4)  is  discontinuous  across  the 

body  surface,  C  being  equal  to  1  outside  the  body  and  to  0  inside.  This 

discontinuity  in  the  value  of  C  on  the  left  side  of  equation  (2.4)  evidently 

is  accompanied  by  a  corresponding  discontinuity  on  the  right  side  of  the 

equation.  Specifically,  the  latter  discontinuity  stems  from  the  integral 

j  k  ifiG^da  representing  the  potential  induced  by  a  surface  distribution  of 

normal  doublets  of  strength  4>  over  b,  as  is  well  known.  An  integral  identity 

valid  for  any  point  £  —  either  outside,  inside  or  exactly  on  the  body 

surface  —  can  be  obtained  by  eliminating  the  discontinuity  in  the  value  of 

C  in  equation  (2.4).  This  can  readily  by  achieved  by  adding  the  term  C.  <J> . 

i 

on  both  the  left  and  right  sides  of  equation  (2.4)  with  given  by 
r  2  f 

C.  =  V  Gdv  =  -  G  da  ,  (2.6) 

1  ->d.  ■'b  n 

i 

as  may  be  obtained  by  using  the  divergence  theorem.  Equation  (2.4)  then  becomes 

<l>*  -  |  GV24.dv  -  |  [G4n-(<|>-<t>*)Gn]da  (2.7) 

since  we  have  C+C^l.  The  integral  identity  (2.7)  is  valid  for  any  point  S 
either  outside,  inside,  or  exactly  on  the  body  surface.  This  new  identity 
thus  is  essentially  equivalent  to  the  set  of  the  three  classical  identities 
(2.5a,b,c),  which  are  valid  exclusively  for  £  outside,  inside,  or  on  the 
body  surface,  respectively.  As  a  matter  of  fact,  these  three  identities 

can  be  obtained  from  the  identity  (2.7)  by  noting  that  we  have  -  /  G  da  = 

(2  -t  b  n 

J ^  V  Gdv  =  0,1,  or  1/2  for  £  outside,  inside,  or  exactly  on  the  body. 

Integral  identities  analogous  to  identities  (2.5a,b,c)  and  2.7)  can 

evidently  also  be  obtained  for  the  "interior  potential-flow  problem". 

These  identities  are  listed  in  appendix  A.  Another  classical  integral 

identity,  obtained  by  combining  these  basic  identities  for  the  exterior  and 

interior  problems,  is  also  given  in  the  appendix. 


ilHiTtrini  ' 
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3.  Basic  integral  equation  of  potential  flow  about  a  body 

2 

In  the  particular  case  of  potential  flow  about  a  body,  we  have  V  <p*0 
in  the  (fluid)  domain  d  outside  the  body,  and  the  integral  identity  (2.7) 
becomes 

?*=  -  |  G*Qda  +  j  (<M*)Gnda  ,  (3.1) 

which  is  an  integral  equation  for  determining  the  velocity  potential  <f>  on  the 
surface  b  of  the  body.  More  precisely,  equation  (3.1)  is  a  nonhomogeneous 
(second  kind)  Fredholm  integral  equation,  of  the  form 


4>(f)  =  f  (t)  +  L(t;$) 


(3.2) 


where  f(£)  is  the  known  (since  <J>  is  specified  on  b)  nonhomogeneous  term  given 
by 


f(t)  -  -  j  G(x,f)<|>n(x)da(x) 


(3.2a) 


and  L(C;<t>)  is  the  linear  transform  of  <{>  defined  by  the  integral 


L(t;<f)  = 


[<j>  (x)  —  <t>  (C  )  ] Gn  (x , \ )  da (x) 


(3.2b) 


The  above  integral  equation  has  previously  been  obtained  by  Chow,  Hou, 
and  Landweber  (1976),  equation  (19),  by  using  a  technique  due  to  Landweber  for 
removing  the  singularity  in  the  dipole  integral  in  the  classical  integral  equation 


=-  G<t>  da  +  <f>G  da  (3.3) 

2  >b  n  Jb  ° 

2 

given  by  the  Green  identity  (2.5c)  in  the  particular  case  when  V  <t>=0  in  d. 

The  manner  in  which  the  integral  equation  (3.2)  was  obtained  in  the  present 
study  supplements  the  Landweber  derivation  and  provides  additional  insight  into 
this  integral  equation.  In  particular,  it  has  been  shown  that  equation  (3.2) 
holds  not  only  for  t  on  the  surface  b  of  the  body,  but  also  for  £  in  the  fluid 
domain  d  outside  the  body.  This  means  that  the  integral  equation  (3.2)  can 


'**»**U 
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in  principle  be  used  to  determine  the  potential  4>  in  the  entire  solution 
domain  d+b,  for  instance  by  using  an  iterative  procedure  based  on  a  recur¬ 
rence  relation  such  as  that  given  by  equation  (3.6).  As  a  matter  of  fact, 
such  an  approach  may  be  necessary  in  the  case  of  subsonic  compressible  flows 
governed  by  a  weakly-nonlinear  equation  of  the  form  V  <j>=g(<j>).  However,  in 
the  case  of  incompressible  flows  it  would  usually  be  much  simpler  to  solve  for 
<j>  on  the  body  surface  b,  and  —  in  the  event  (rare  in  reality)  that  knowledge 
of  <(>  outside  b  is  in  fact  required  —  to  determine  4>  outside  b  by  means  of 
equation  (2.5a),  which  here  takes  the  simplified  form 


' ) 


G<j)  da 

n 


<j>G  da 
n 


(3.4) 


Equation  (3.2)  also  holds  for  £  inside  the  body.  It  thus  might  appear  that 
this  integral  equation  must  also  define  the  potential  $  inside  the  body.  This 
result,  were  it  true,  would  certainly  be  quite  surprising,  indeed  fundamentally 
unacceptable,  for  it  would  mean  that  the  "exterior  boundary-value  problem" 
defined  by  the  Laplace  equation  in  the  exterior  domain  d  and  a  Neumann  condition 
on  b  would  define  a  solution  in  the  interior  domain  d^.  It  can  easily  be  shown 
however,  that  equation  (3.2)  allows  the  potential  4>  to  be  extended  inside  b 
in  an  entirely  arbitrary  manner.  Indeed,  equation  (3.1)  can  be  written  in 
the  form 


rj 

i 

to 


C4> 


* 


G(J>  da 
n 


(3.5) 


where  C  is  given  by  C=l+/ljGnda.  We  have  CsO  for  £  inside  b,  as  may  easily  be 
verified  (and  indeed  has  already  been  shown  in  the  previous  section) ,  so  that 
equation  (3.5)  clearly  does  not  define  inside  b.  For  £  outside  b,  on  the 
other  hand,  we  have  C=l,  and  equation  (3.5)  becomes  equation  (3.4),  which 
clearly  defines 

An  interesting  feature  distinguishing  the  Landweber  integral  equation 
(3.2)  from  the  classical  integral  equation  (3.3)  is  that  whereas  the  dipole 
integral  in  the  usual  integral  equation  (3.3)  is  discontinuous  at  the  body 
surface  b  (as  was  discussed  in  the  previous  section,  and  is  indeed  well  known), 
the  corresponding  integral  in  equation  (3.2)  is  a  continuous  function  (the 
dipole  strength  vanishes  as  the  point  of  integration  x  and  the  field 


point  £  coincide).  As  a  matter  of  fact,  the  factor  1/2  on  the  left  side  of 

equation  (3.3)  is  correct  only  at  points  £  where  b  is  smooth  (that  is,  has 

a  tangent  plane),  as  was  also  discussed  in  the  previous  section.  The  classical 

integral  equation  (3.3)  thus  requires  evaluation  of  a  discontinuous  function 

exactly  on  the  surface  of  discontinuity  of  that  function.  This  awkward  problem, 

notably  from  the  point  of  view  of  numerical  calculations,  is  avoided  in  the 

integral  equation  (3.2).  As  a  matter  of  fact,  the  integrand  (d>  — ^.) G  in  the 
->■  n 
integral  L(£;$)  defined  by  equation  (3.2b)  is  non-singular  (i.e.  remains  finite) 

as  x  -*•  £  (as  can  easily  be  verified,  and  indeed  is  shown  in  Chow,  Hou,  and 

Landweber) ,  which  evidently  is  advantageous  for  purposes  of  numerical  calculations; 

for  two-dimensional  flows,  the  integrand  (iHOG  actually  vanishes  as  x  -*■  £. 

A  choice  of  methods  is  available  for  solving  the  nonhomogeneous  Fredholm 

integral  equation  (3.2).  In  particular,  a  natural  method  of  solution  consists 

in  using  an  iterative  procedure.  An  obvious  recurrence  relation  is  that  obtained 

(k)  (k+1) 

by  simply  replacing  <j>  by  <f>  and  </>  on  the  right  and  left  sides,  respectively, 

of  equation  (3.2),  that  is 

<l>(k+1)(£)  =  f (£)  +  L(£;$(k)),  k>0  ,  (3.6) 


where  the  initial  (zeroth)  approximation  must  be  specified  somehow.  This 

recurrence  relation  has  in  fact  been  used  by  Chow,  Hou,  and  Landweber  (1976), 

equation  (20),  with  the  initial  approximation  4>  ^  (£)  taken  as  twice  the  given 

term  f(£).  The  practical  usefulness  of  such  an  iterative  method  of  solution 

(k) 

depends  on  fast  convergence  of  the  successive  iterative  approximations  <j>  , 

which  in  turn  depends  on  selection  of  a  reasonably-good  initial  approximation 
<f>^\  Various  initial  approximations  and  associated  iterative  approximations 
will  be  investigated  in  this  study. 


y.t-turN'i 


-10- 


4.  Potential  flow  about  ellipses  and  ellipsoids 

In  this  section,  the  integral  equation  (3.2)  and  the  associated  recurrence 

relation  (3.6)  are  investigated  by  considering  the  simple  cases  of  potential 

flows  due  to  translatory  motions  of  elliptical  cylinders  and  ellipsoids.  Two- 

dimensional  potential  flows  due  to  translatory  motions  of  an  ellipse  in  the 

directions  of  its  major  and  minor  axes  are  considered  first.  The  ellipse  is 

2  2  2 

defined  by  the  equation  x  +y  /b  -1,  where  b  is  the  thickness  ratio  of  the 
ellipse  (b*thickness/chord) .  Both  "longitudinal"  and  "transversal"  motions  of 
the  ellipse  with  unit  velocity  in  the  x  and  y  directions  of  the  major  and  minor 
axes,  respectively,  are  considered. 

For  longitudinal  motion,  the  velocity  potential  on  the  surface  of  the 
ellipse  is  given  by 

$  =  -bx  .  (4.1) 

Furthermore,  it  can  be  shown  that  the  given  term  f  defined  by  equation  (3.2a) 
actually  is  proportional  to  the  potential  ;  specifically,  we  have 

f/4>  *  1/ (1+b)  .  (4.1a) 

Use  of  equation  (4.1a)  in  the  integral  equation  (3.2)  then  yields 

L (<f>)/<|>  -  b/(l+b) .  (4.1b) 

For  transversal  motion,  on  the  other  hand,  we  have 

♦  =  -(l-x2)1/2  (4.2) 

on  the  upper  portion  (y>0)  of  the  surface  of  the  ellipse.  The  terms  f  and  L($) 
in  the  integral  equation  (3.2)  are  also  proportional  to  specifically,  we 
have 


f/$  -  b/(l+b)  ,  L(*)/$  -  1/ (1+b)  .  (4.2a,b) 

It  thus  can  be  seen  that  the  value  of  f/$  for  longitudinal  (transversal)  motion 
of  the  ellipse  is  identical  to  the  value  of  L(4)/4>  for  transversal  (longitudinal) 
motion. 

The  functions  f/$  and  L($)/<|>  for  longitudinal  (L)  and  transversal  (T) 
motions  are  represented  in  figure  2  for  values  of  b  varying  between  0  and  1, 
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for  which  the  ellipse  becomes  a  flat  plate  and  a  circle,  respectively.  Figure  2 
clearly  shows  that  the  known  term  f  in  the  integral  equation  (3.2)  is  larger 
than  the  unknown  term  L(<j>)  for  longitudinal  motion  of  the  ellipse,  whereas  the 
opposite  is  true  for  transversal  motion;  in  the  intermediate  case  of  a  circle 
(b=l),  we  have  f *L ($ )  =<J> /2 .  In  the  particular  case  of  a  thin  ellipse  (b<<l), 
equations  (4.1)  show  that  we  have 

<(>  =  0(b),  f  =  0(b),  L(4>)  =  0(b2)  and  L(4>)  «  f  =  )  (4.3a) 

for  longitudinal  motion,  whereas  for  transversal  motion  equations  (4.2)  yield 

<f>  =  0(1),  f  =  0(b),  L(<J>)  «  0(1)  and  f  <<  L(<j>)  *  <fr  .  (4.3b) 

The  given  term  f  in  the  integral  equation  (3.2)  may  thus  be  used  as  a  "thin- 
wing  approximation"  for  longitudinal  motion  of  the  ellipse,  but  clearly  not  for 
transversal  motion. 

Three-dimensional  potential  flows  due  to  translatory  motions  of  an  ellipsoid 

2  2  2  2  2  2 

are  now  considered.  The  ellipsoid  is  defined  by  the  equation  x  +y  /6  +z  /B  y  =1 , 

where  the  parameters  6  and  y  may  take  values  between  0  and  1,  that  is  0<6<1  and 

0<_y<1.  Translatory  motions  of  the  ellipsoid  with  unit  velocity  in  the  x,y,  and 

z  directions  of  the  major,  intermediate,  and  minor  axes  are  considered.  The 

velocity  potentials  corresponding  to  these  three  fundamental  translatory  motions 
x  y  z 

are  denoted  by  $  ,4>J ,  and  ij>  ,  and  the  corresponding  potentials  f  given  by 
equation  (3.2a)  are  similarly  denoted  by  fX,f^,  and  fZ. 

For  translation  with  unit  velocity  in  the  positive  x  direction,  we  have 
4>n=i*n,  where  i  is  the  unit  vector  along  the  positive  x  axis,  and  n  in  the  unit 
inward  normal  vector  to  the  surface  (b)  of  the  ellipsoid  as  previously  defined. 
Equations  (3.2a)  and  (2.1)  then  yield 

^  ^  fb  * 
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where  is  the  domain  inside  the  surface  of  the  ellipsoid.  By  defining  the 
function  F(£)  as 

f<?>  "  iir  I  -rrr  dv(S)  5  I  o<x,t) dv<5)  ,  (4.4) 

4’  |J-!l  JaI 

x 

we  may  finally  express  the  potential  f  in  the  form 

f*(t)  =  3F(t)/3C  .  (4.5a) 


Similarly,  we  have 

fy(!)  -  3F(f)/3n  ,  f Z(t)  -  3F(t)/3S  .  (4.5b,c) 

X  v 

It  is  known  from  Havelock  (1931)  that  the  velocity  potentials  <J>  ,4>  ,  and 
<5>  for  unit-velocity  translations  of  the  ellipsoid  along  the  x,y,  and  z  axes 
can  be  expressed  in  terms  of  the  function  F(t)  defined  by  equation  (4.4). 
Specifically,  we  have 

d>X  -  (l/AX)3F/3£,  $y  =  (l/Ay)3F/3n,  <t>Z  =  (1/AZ)3F/3C  ,  (4.6a,b,c) 

x  y  z 

where  A  ,A  ,  and  A  are  constants  defined  by  the  integrals 


g  ydt 
2U+1)E 


g2ydt 

2(t+S2)E 


2(t+62y2)E 


(4. 7a,b,c) 


2  22  1/2 

in  which  E  is  the  expression  given  by  E»[ (t+1) (t+6  ) (t+6  y  )] 

It  may  then  be  seen  from  equations  (4.5)  and  (4.6)  that  the  given  term 
f(C)  in  the  integral  equation  (3.2)  is  proportional  to  the  velocity  potential 
<J>(£)  for  translatory  motions  of  an  ellipsoid.  Specifically,  we  have 


AX,  fy/*y  -  Ay,  fZ/*Z  -  A2 


(4.8a,b,c) 


where  the  constants  Ax,Ay,  and  Az  are  defined  by  the  integrals  (4.7).  These 
constants  of  proportionality  are  represented  in  figure  3  as  functions  of  the 
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parameters  0  and  y  defining  the  ellipsoid.  This  figure  is  subdivided  into 
four  parts. 

The  lower  left  part  of  figure  3  represents  AX  as  a  function  of  8(0<B<1) 
for  given  values  of  y(y=0, . 2, .4, .6, . 8,1) .  We  have  l>Ax>2/3,  wit1'  AX=2/3  in 
the  limiting  case  of  a  sphere  (0=1=y) •  For  translation  of  an  ellipsoid  along 
its  major  axis,  the  known  term  f  in  the  integral  equation  (3.2)  may  then  be 
seen  to  be  the  dominant  term  in  comparison  with  the  unknown  term  t,(4>),  which 
does  not  exceed  <f>/3.  In  particular,  the  value  of  the  constant  Ax  is  quite 
close  to  1  for  small  values  of  0.  This  shows  that  the  term  f  provides  a  fairly 
good  "slender-body  approximation"  to  the  velocity  potential  <J>  for  translation 
of  a  slender  ellipsoid  in  the  direction  of  its  major  axis.  Comparison  between 
the  values  of  A=f/<f>  and  of  A  =f  /<£  represented  in  figures  2  and  3  for  longi¬ 
tudinal  motions  of  an  ellipse  and  of  an  ellipsoid,  respectively,  show  that  we 

X  X 

have  A  (b) < A  (8;y=1)£A  (8;0<y^1)  for  corresponding  values  of  the  thickness  and 
slenderness  parameters  b  and  0,  respectively.  More  precisely,  we  have 
l/2<A(b)=l/(l+b)<l,  and 


-|  <_  AX(0;y=1)  = 


1-8 


[1+ 


2 (1-02) 1/2 


l+(l-&  )  ' 


as  may  be  obtained  by  evaluating  the  integral  (4.7a)  in  the  special  case  of  a 
spheroid  (y~1) •  We  thus  have,  for  instance,  A=.833  for  b=.2  and  .944£AX<  1  for 
6=  .2  and  1^Y^_0 . 

The  lower  right  and  upper  left  parts  of  figure  3  correspond  to  translation 

of  the  ellipsoid  in  the  direction  of  its  intermediate  axis,  and  represent  the 

constant  Ay(B,Y)=fy/4>y  both  as  a  function  of  8  for  given  values  of  y  (lower 

right  part  of  figure  3)  and  as  a  function  of  y  for  given  values  of  0  (upper 

left  part).  It  may  be  seen  that  we  have  l>Ay^l/2,  with  Ay=l/2  in  the  limiting 

case  0=0, y*1  corresponding  to  a  circular  cylinder.  Finally,  the  upper  right 

part  of  figure  3  corresponds  to  translation  of  the  ellipsoid  along  its  minor 

z  z  z 

axis,  and  represents  A  (8,y)“  f  / d>  as  a  function  of  y  for  given  values  of  0. 

In  this  case  we  have  0^AZ<^  2/3,  with  Az=2/3  for  B“1”Y  (sphere).  In  summary, 
the  given  term  f  in  the  integral  equation  (3.2)  may  be  seen  to  be  dominant  for 
translation  of  the  ellipsoid  in  the  directions  of  its  major  and  minor  axes, 
whereas  for  translation  along  its  minor  axis  the  term  L(<f>)  is  usually  dominant, 
except  for  values  of  y  approximately  given  by  1-8/2<y.5  1  • 


I 


We  conclude  this  investigation  of  the  integral  equation  (3.2)  for  flows 
due  to  translatory  motions  of  an  ellipsoid  by  examining  the  iterative  procedure 
defined  by  the  recurrence  relation  (3.6)  and  the  initial  approximation  <(>^=0. 

By  using  the  fact  that  we  have  f-A$,  we  may  obtain  <J>^^=f  ,<(/^“[l+(l-X)  ]f 
[1+(1-A)+(1-A) 2]f  t  and  more  generally  (J>^n^=[l+(1-A)+. .  .  +  (1-A)n  ^]f  for  n>l.  This 
series  converges  to  f /  [  1—  ( 1—  A)  ]  =  f /X =4>  as  n-*»  provided  -1<1-A<1.  We  have  1^A>0, 
which  yields  0£l-A<:l.  The  above-defined  series  and  associated  iterative  procedure 
therefore  converge,  except  in  the  limiting  case  A=0  corresponding  to  translation 
of  an  elliptical  disk  (y=0)  in  the  direction  normal  to  its  plane  (direction  of 
the  minor  axis  of  the  degenerate  ellipsoid). 

More  precisely,  the  relative  error,  say,  associated  with  the  n^-  itera¬ 
tive  approximation  defined  as  e  =  (4>— 4>  ,  can  be  shown  to  be  given  by 

n  n 

en=(l-A)  .  The  number  of  iterations  n  required  for  obtaining  an  approximation 
to  the  potential  within  a  specified  relative  error  e  then  is  given  by  n>£ne/£n(l-A) . 


For 

instance. 

for  e*10  ^  the 

following  number  of 

iterations 

is  required 

A 

.9 

.8 

.7 

.6  .5 

.4 

.3 

.2 

.1 

n 

2 

3 

4 

5  7 

9 

13 

21 

44 

depending  on  the  value  of  A.  As  an  application,  we  consider  translatory  maneuver¬ 
ing  motions  of  an  ellipsoidal  ship  hull  form  with  values  of  f?  and  y  equal  to  .13 
and  .7,  corresponding  to  the  fairly  typical  values  of  the  beam/length  and  draft/ 
length  ratios  .15  and  .05,  respectively. Figure  3  indicates  corresponding  values 
of  A  *f  /tf>  and  A  -f7/(jr  approximately  equal  to  .975  and  .6,  for  which  we  need 
2  and  5  iterations,  respectively. 


i 
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5.  Potential  flow  about  ogives 

The  remarkable  property,  enjoyed  by  ellipsoids,  that  the  given  term 
f  in  the  integral  equation  (3.2)  is  proportional  to  the  velocity  potential 
<t>  does  not  hold  in  general  of  course.  It  thus  is  useful  to  consider  other 
body  shapes  besides  ellipsoids.  For  simplicity,  calculations  have  been  per¬ 
formed  for  two-dimensional  potential  flows  due  to  translations  of  ogives  along 
their  major  and  minor  axes.  A  main  reason  for  considering  flows  due  to 
translations  of  ogives  is  that  the  velocity  potential  $  can  be  determined 
analytically,  by  using  conformal  mapping.  Another  recommendation  of  ogives 
is  that  their  pointed  shapes  obviously  supplement  ellipses. 

The  ogives  are  defined  by  the  equation 


x2  +  [|y|+(l-b2)/2b]2  =  [(l+b2)/2b]i 


(5.1) 


where  b  is  the  thickness  ratio  (0<b<l)  and  -l<x<l,or  by  the  parametric  equations 


x  = 


l-t‘ 


l+t2+2t(l-b2)/(l+b2) 


4bt/ (1-H>  ) 


l+t2+2t(l~b2)/(l+b2) 


(5.1a,b) 


\  jp  2  —1 

where  t  is  defined  as  t=(tanj)  with  6=2(1-—  tan  b) ,  and  the  parameter  X 
varies  between  the  values  0  and  tt/2  (we  have  |x|=l  and  y=0  for  X=0,  and  x=0 

and  |y|*b  for  X=tt/2).  In  the  limit  b=l,  the  ogive  becomes  the  unit  circle 

2  2.  2 
x  +y  =1,  while  in  the  limit  b-*-0  we  have  |y|M>(l-x  ). 

The  velocity  potentials,  on  the  body  surface  b,  due  to  unit-velocity 

2  1/2 

translations  of  ellipses,  defined  by  the  equation  |y[=b(l-x  )  ,  and  ogives 

with  corresponding  values  of  the  thickness  ratio  b  are  represented  in  figures 
4  and  5  for  longitudinal  and  transversal  motions  (translations  along  the  major 
and  minor  axes),  respectively.  These  figures  represent  the  potential  4>  for 
x  varying  between  the  forward  stagnation  point  and  the  midchord  point,  that 
is  for  l>_x>0  and  0<x<l  in  figures  4  and  5,  respectively.  The  potential  <j>  on 
the  surface  of  an  ellipse  in  longitudinal  motion  is  a  linear  function  of  x; 
more  precisely  we  have  $=-bx,  as  is  given  by  equation  (4.1).  The  potential 
on  the  surface  of  an  ogive  is  given  by 


x  -(cosX)/(l-—  tan  *b) 

TT 


(5.2a) 


where  x(0£x<l)  is  given  by  equation  (5. la), and  0£X<tt/2.  This  potential  is 


-18- 


no  t  a  monotonic  function  of  x,  as  may  be  seen  from  figure  4.  The  potential 

<f>  on  the  upper  surface  (y>0)  of  an  ellipse  in  transversal  motion  is  given  by 
2  1/2 

<|>*-(l-x  )  ,  for  any  value  of  the  thickness  ratio  b.  For  transversal  motion 

of  an  ogive,  we  have 

<t>  =  y  -(sinX)/(l-  ~  tan”H>)  ,  (5.2b) 

where  y(b>y>0)  is  given  by  equation  (5.1b),  and  0<X<tt/2.  This  potential  be- 

2  1/2 

comes  identical  to  the  potential  <l=-(l-x  )  '  for  an  ellipse  in  the  limiting 
cases  b=l  and  b*0  (corresponding  to  flow  past  a  circle  and  a  normal  flat 
plate, respectively ,  as  is  indicated  in  figure  5. 


E 


Figure  5  -  The  velocity  potential  $  on  the  surfaces  of 

ellipses  and  ogives  in  transversal  translation 


*■  i*war\ 
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6.  Simple  explicit  approximations  for  the  potential 

It  was  shown  in  section  4,  equations  (4.1a)  and  (4.2a),  that  the  velocity 
potential  $  on  the  surface  of  an  ellipse,  with  thickness  ratio  b,  is  related  to 
the  given  term  f  in  the  integral  equation  (3.2)  by  the  equations 

<t>  =  (1+b)  f  and  <J>  =  (1+1/b)  f  (6.1a,b) 

for  longitudinal  and  transversal  translations,  respectively.  This  result  suggests 
that  useful  simple  explicit  approximations  to  the  velocity  potential  for  longitudinal 
and  transversal  translations  of  other  cylinders  might  be  given  by  the  potentials 
(l+b)f  and  (1+1/b) f,  respectively,  where  b  is  the  thickness  ratio  of  the  cylinder 
and  f  is  defined  by  formula  (3.2a). 

In  particular,  these  explicit  approximations  might  be  used  as  zeroth-order 
(initial)  approximation  <f>  ^  in  the  recurrence  relation  (3.6),  that  is 

=  (l+b)f  and  $(0)  =  (1+1/b)  f  (6. 2a,b) 

for  longitudinal  and  transversal  motions,  respectively.  The  first-order  iterative 
approximations  given  by  equations  (3.6)  and  (6.2)  are 

=  f  +  (l+b)L(f)  and  <i(1)  =  f  +  (l+l/b)L(f)  (6.3a,b) 

for  longitudinal  and  transversal  motions,  respectively.  These  approximations  are 
also  exact  for  translations  of  elliptical  cylinders,  as  may  be  verified  by  using 
equations  (6.1a,b)  in  the  term  L (<p )  on  the  right  side  of  the  integral  equation 
(3.2). 

The  approximations  (6.2b)  and  (6.3b)  for  transversal  translation  lead  to 

numerical  difficulties  in  the  limiting  case  of  thin  cylinders,  that  is  for  b<<l. 

Indeed,  we  have  f-K)  as  b-v0,  so  that  the  terms  f/b  and  L(f)/b  are  in  the  indeterminate 

form  0/0.  An  alternative  to  these  approximations,  for  the  case  of  transversal 

translation  of  thin  cylinders,  consists  in  using  the  potential,  say,  for  flow 

(0)  Ip 

past  a  normal  flat  plate  as  zeroth-order  approximation  $  .  This  then  yields 

the  zeroth-and  first-order  approximations 

=  iifp  *  -(l-x2)1/2  for  y  >  0,  and  =  f  +  L(<f> f p )  - 


(6.4a,b) 
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These  approximations  again  are  exact  in  the  particular  case  of  elliptical  cylinders, 
as  may  readily  be  seen  from  equations  (4.2)  and  (3.2). 

The  above-defined  approximations  are  compared  to  the  exact  potential  «J>  on 
the  surfaces  of  ogives  in  longitudinal  and  transversal  translations  in  figures  6, 

7,  and  8.  Figure  6  shows  the  potential  <p  and  the  approximations  f,  (l+b)f,  and 
f+(l+b)L(f)  for  longitudinal  motion  of  thin  ogives,  with  thickness  ratio  b  equal 
to  .1,  .2  and  .3.  This  figure  clearly  shows  that  multiplication  of  the  basic  thin- 
body  approximation  f  by  the  constant  (1+b)  yields  significant  improvements,  and 
the  improvements  obviously  are  more  important  for  larger  values  of  b.  Figure  6 
also  depicts  the  improvements  of  the  zeroth-order  approximation  <j>^^  =  (l+b)f  due 
to  the  first-order  approximation  <j>  ^  =f+(l+b)  L(f )  .  These  improvements  are  apparent 
also  from  figure  7,  in  which  the  approximations  (l+b)f  and  f+(l+b)L(f)  are  compared 
to  the  exact  potential  4>  for  longitudinal  translation  cf  ogives  with  thickness 
ratio  b  equal  to  .25,  .5,  .75,  and  1.  For  b=l,  the  ogive  is  a  circle  and  the  approxi¬ 
mations  (l+b)f  and  f+(l+b)L(f)  are  identical  to  ()>5-x. 

Figure  8  corresponds  to  transversal  translation  of  ogives  and  represents 

the  potential  <p  and  the  approximations  (l+l/b)f,  f+(l+l/b)L(f )  and  f+L((J>f  )  for 

2  1  7  0 

b=0,  .1,  .25,  .5,  .75,  and  1;  the  zeroth-order  approximation  <{>£  =-(l-x  )'  L  has 

already  been  represented,  together  with  the  potential  ij> ,  in  figure  5.  In  the 

limiting  case  b=l,  all  these  approximations  are  identical  to  the  exact  potential 

2  1/2 

<P,  which  is  given  by  <f)=—  ( 1— x  )  .  For  fairly  large  values  of  b,  say  for  b=.75 

and  for  b=.5,  the  zeroth-and  first-order  approximations  (l+l/b)f  and  f+(l+l/b)L(f) 
may  be  seen  to  be  fairly  close  to  the  potential  (j> ,  and  somewhat  superior  to  the 
alternative  approximations  <|>fp  and  f+L(<t>fp).  Figure  8  also  shows  that  the  first 
approximation  f+(l+l/b)L(f)  improves  upon  the  zeroth  approximation  (l+l/b)f.  The 
approximations  (l+l/b)f  and  f+(l+l/b)L(f )  may  however  be  seen  to  deteriorate  as  b 
decreases,  and  the  alternative  approximations  and  f+L^^)  clearly  are  superior 
for  fairly  small  values  of  b,  say  for  b<.25.  In  particular,  in  the  limiting  case 
of  a  flat  plate  (b=0) ,  the  approximations  <j>,  and  f+L($,  )HL(<{>,  )  become  exact. 


2.1/2 


fp'^^fp' 


equal  to  the  potential  <j>=-(l-x  )  .  On  the  other  hand,  we  have  (1+1/b) fvf /b 

and  f+(l+l/b)L(f ),vL(f /b),'.f /b  as  b+0,  so  that  the  approximations  (1+1/b) f  and 
f+(l+l/b)L(f)  become  identical  in  the  limit  b=0;  however,  the  limit  of  f/b  as 
b>0  is  not  equal  to  the  potential  ,  as  may  be  seen  from  figure  8. 


re  8  -  The  potential  $  and  the  approximations  (l+l/b)f,  f+(l+l/b)L(f ) , 

and  f+L((J>,  )  on  the  surfaces  of  ogives  in  transversal  translation 
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7.  A  modified  integral  equation  and  related  explicit  approximation 


The  preceding  investigation  of  flows  due  to  translations  of  ellipsoids  and 
ogives  suggests  that  it  may  be  advantageous  to  express  the  solution  $(£)  of  the 
integral  equation  (3.2)  in  the  form  4> (£)=k(£) f (£) ,  where  f(£)  is  the  potential 
given  by  formula  (3.2a)  and  k (£)  is  the  function  defined  as  k=<J>/f.  Figures  6, 

7,  and  8  for  flows  due  to  translations  of  ogives  suggest  that  k(£)  may  generally 
be  expected  to  be  a  slowly  varying  function.  This  function  in  fact  is  a  constant 
in  the  particular  cases  of  translations  of  ellipsoids,  as  was  shown  previously 
and  is  specifically  indicated  by  equations  (4.8).  More  generally,  we  will  express 
the  potential  <j>0t)  in  the  form  <j>  (t)=k(CH(£)  ,  where  is  some  supposedly- 

given  potential.  An  obvious  choice  for  sfi  is  the  given  term  f  in  the  integral 
equation  (3.2),  as  was  just  discussed.  Other  choices  for  the  potential  ^  can 
however  be  made.  For  instance,  for  two-dimensional  flow  due  to  transversal  trans¬ 
lation  of  a  thin  cylinder,  the  potential  <p  could  be  taken  as  the  potential 
2  1/2 

<j>£p=-(l-x  )  for  flow  past  a  normal  flat  plate,  as  was  discussed  in  the 
previous  section. 

By  expressing  the  term  4>  (x)-<ji(£)=k(x)iJj(x)-k(0'K£)  in  equation  (3.2b)  in 
the  form  k(f)  [<Hx)-iK?)  ]+[k(x)-k(£)  ]iJj(x)  ,  the  integral  equation  (3.2)  becomes 


[k(x)-k(£) ]^(x)Gnda. 


<i>(£)  =  f(£)  +  k(£)  [iKx)->K£)]G  da  + 

'b 


Multiplication  of  this  equation  by  ifi(£)  and  use  of  the  relation  kiji=<}>  then  yields 
<f>(t)[<Kt)  -  |  [<Kx)-Mt)]Gnda]  =  f (t)<K£)  +|  [<j>(x)iKf)-lj>(£)4'(x))Gnda  . 


1  b  “  J  '  b 

We  thus  may  obtain  the  modified  integral  equation 

<f>(t)  =  </>(£)  +  L  ’(£;4>)  , 


(7.1) 


where  the  potential  $(£)  is  defined  as 


<fKO  -  f(£)^(£)/[iK£)-L(£;^)]  , 


(7.1a) 


with  L(£;iJ))  *  /.  [^(x)->(j(£)  ]G  (x,£)da(x)  as  is  given  by  equation  (3.2b),  and  the 
d  n 

term  L’(£;$)  is  the  linear  transform  of  $  defined  by  the  integral 


L'(£;4>)  -  j  [<i(x)^(£)-<()(t)tKx)  ]Gnday,[i^(t)-L(£  ;\p)  ]  . 


(7.1b) 
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The  term  L*  ;4>)  obviously  vanishes  if  the  function  ip<0  is  proportional  to  the 
potential  <t>(t),  whereas  the  function  $(t)  then  becomes  identical  to  $(t),  as  may 
be  verified  by  using  the  integral  equation  (3.2)  in  expression  (7.1a). 

It  then  follows  from  equations  (4.8)  that  the  potential  defined  by  equation 
(7.1a)  with  fp  taken  as  the  term  f  given  by  formula  (3.2a),  that  is  the  potential 
ij)  given  by 

<P(t)  =  f2(l)/ff(|)-L(|;f)]  ,  (7.2) 

is  exact  in  the  particular  cases  of  potential  flows  due  to  translations  of  ellipsoids. 
In  general,  formula  (7.2)  defines  an  explicit  approximation  for  the  velocity 
potential  <t> .  This  approximation  is  compared  to  the  exact  potential  tp  on  the 
surfaces  of  ogives  in  longitudinal  and  transversal  translations  in  figures  9,  10, 
and  11. 

2 

Figure  9  shows  the  potential  $  and  the  approximations  f,f+L(f),  and  f  / [ f — L ( f ) ] 

for  longitudinal  motion  of  thin  ogives,  with  thickness  ratio  b=  .1,  .2,  and  .3. 

The  approximations  f  and  f+L(f)  are  the  first  and  second  iterative  approximations 

associated  with  the  straightforward  recurrence  relation  (3.6)  and  the  initial 

approximation  Computationally,  the  approximations  f+L(f)  and  f2/[f-L(f)] 

2 

are  equivalent.  However,  the  approximation  f  /[f-L(f)]  may  be  seen  to  be  superior 

2 

to  the  approximation  f+L(f).  The  superiority  of  the  approximation  f  / [ f — L ( f ) ] 
in  comparison  with  the  approximation  f+L(f)  is  apparent  also  from  figure  10, 
where  these  two  approximations  are  compared  to  the  exact  potential  <p  for  longi¬ 
tudinal  translation  of  ogives  with  thickness  ratio  b=.25,  .5,  .75,  and  1.  In 

2 

particular,  the  approximation  f  / [ f— L( f ) ]  is  identical  to  the  potential  $=-x 

in  the  limiting  case  b=l,  corresponding  to  a  circular  cylinder;  for  b=.75,  this 

approximation  is  practically  indistinguishable  from  4>  on  the  scale  of  figure  10. 

2 

It  is  also  interesting  to  compare  the  approximation  f  / [ f — L ( f ) ]  to  the 
(computationally-equivalent)  approximation  f+(l+b)L(f) •  Comparison  of  figures 
10  and  9  to  figures  7  and  6,  respectively,  show  that  the  approximation  f  / [ f — L( f ) ] 
is  superior  to  the  approximation  f+(l+b)L(f)  except  for  sufficiently-thin  ogives, 
say  for  b<^2,  for  which  the  two  approximations  are  comparable.  Finally,  figures 
10,9,7,  and  6  incidentally  provide  a  comparison  between  the  straightforward 
second  approximation  f+L(f)  and  the  modified  first  approximation  (l+b)f. 

Figures  6  and  9  show  that  the  approximation  f+L(f)  is  somewhat  superior  to  the 
approximation  (l+b)f  for  b=.l,  and  that  these  two  approximations  are  comparable 
for  b*.2  and  .3.  However,  for  thicker  ogives,  say  for  . 5<b<l,  figures  7  and  10 
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show  that  the  (computationally-simpler)  approximation  (l+b)f  in  fact  is  superior 
to  the  approximation  f+L(f). 

Figure  11  corresponds  to  transversal  translation  of  ogives  with  thick¬ 
ness  ratio  b=0,  .1,  .25,  .5,  .75,  and  1.  This  figure  shows  the  approximation 

fi{>,  /[<j>f  — L (4> f  )]»  obtained  by  selecting  the  potential  if  in  equation  (7.1a)  as 
rP  rP  IP  2  1/2 

the  potential  d> £p=— ( 1— x  )  of  flow  past  a  normal  flat  plate,  and  the  computa- 

tionally-equivalent  approximation  f+L(d> ^p) »  obtained  previously  in  equation  (6.4b). 

In  the  limiting  cases  b=0  and  b=l,  both  of  these  approximations  are  identical  to 

the  exact  potential  4>=— (1— x  )  'L  .  More  precisely,  we  have  f=<j> /2  ,  4>  ,L($^  )*4>/2 , 

and  if>£p-L(<J> )  =<f> /2  for  b=l,  while  as  b+0  we  have  f->-0,  and 

as  may  be  verified  from  the  integral  equation  (3.2);  the  approximation 

f'f’fp/[<t,£p-L(if)  jp)  ]  thus  is  in  the  indeterminate  form  0/0  for  b=0,  although  we 

have  f/t<J’fp-L(4>£p)  1+1  as  b+0.  Disregarding  this  potential  numerical  difficulty 

for  very  thin  cylinders,  figure  11  shows  that  the  approximation  fif>^  / [d> jrp— L. ) ] 

is  superior  to  the  approximation  f+L(<{>f  )  • 

Ip  2 

Figure  11  also  shows  the  approximation  f  /[f-L(f)]  for  b=l,  .75,  and  .5. 

In  the  limiting  case  b=l,  we  have  f =<f> /2  and  L(f)=>if>/4,  so  that  we  have  f-L(f)=$/4 

2 

and  the  approximation  f  / [ f — L( f )  ]  is  identical  to  4> .  For  b=.75  and  .5,  this 

approximation  may  be  seen  to  be  fairly  good,  and  roughly  comparable  to  the  approxi- 

mation  f^^/ [<f>fp-L(<f>  fp)  ] .  However,  the  approximation  f  /  [ f — L ( f )  ]  deteriorates  as 

b  decreases,  and  it  is  not  shown  for  b<.5.  The  reason  for  the  fact  that  the 
2 

approximation  f  / [ f — L( f ) ]  is  not  useful  for  transversal  translation  of  thin 

cylinders  is  that  the  terms  f  and  L(f)  then  are  comparable,  as  can  easily  be  shown 

by  considering  the  particular  cases  of  ellipsoids.  For  flows  due  to  translations 

of  ellipsoids  the  given  term  f  in  the  integral  equation  (3.2)  is  proportional  to 

the  exact  potential  <f>,  that  is  we  have  f»A<f>,  as  is  indicated  in  equations  (4.8). 

Use  of  this  relation  in  th~  integral  equation  (3.2)  then  yields  L(  f )  =AL(<}>)=A  (<>-f  )  = 

2  2  2  2  2 

A(l-A)4>,  so  that  we  have  f-L(f)=A  <f>  and  f  / ( f  — L ( f )  1  “A  !>  /\  .  The  approximation 

2 

f  / [ f— L( f ) ]  thus  is  exact,  as  was  indeed  noted  previously,  but  numerical  diffi¬ 
culties  may  clearly  be  expected  if  L(f)=A  (1— X ) <J>  is  comparable  to  f =X ,  which 
occurs  if  A  is  small.  For  practical  purposes,  it  may  be  preferable  not  to  use 
the  approximation  f  /[f-L(f)]  for  values  of  A  less  than  about  .4,  for  which  wo 

have  f*.4$,  L(f)*.24<{> ,  and  f — L (f)«.164>.  Figure  3  then  indicates  that  the  approxi- 
2 

mation  f  / [ f— L( f ) ]  might  be  used  for  translations  of  a  body  along  its  major  and 
intermediate  axes,  as  well  as  along  its  minor  axis  if  the  body  Is  not  too  flat. 

More  precisely,  figure  3  shows  that  we  have  A^.4  for  translation  of  an  ellipsoid 
alongs  its  minor  axis  if  y<.65-.  3  approximately. 
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Figure  10  -  The  potential  and  the  approximations 


f+L(f)  and  f2/[f-L(f)] 


on  the  surfaces  of  ogives  in  longitudinal  translation 
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Figure  11  -  The  potential  4>  and  the  approximations  f+L^^), 

and  f  /[f-L(f)]  on  the  surfaces  of  ogives  in  transversal  translation 
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8.  Conclusion 

The  foregoing  study  of  flows  due  to  translation  of  ellipsoids  and  ogives 
shows  that  simple  explicit  formulas  for  the  velocity  potential  can  provide 
realistic  approximations,  which  may  actually  be  sufficient  for  many  practical 
purposes.  Two  simple  explicit  approximations  appear  to  be  of  particular  interest. 
These  are  the  approximation  f(£)  given  by 


G(x,t)<t>n(x)da(x) 


(8.1) 


as  is  defined  by  formula  (3.2a),  and  the  approximation  $(5)  given  by  formulas 
(7.2)  and  (3.2b),  that  is 


!(t)/[f(t)- 


{f(x)-f (f)}Gn(x,f)da(x) ]. 


(8.2) 


The  approximation  f(£)  is  a  (first-order)  slender-body  approximation  which  may  be 
useful  for  longitudinal  translation  of  a  slender  body,  such  as  a  ship  form  for 
instance.  As  a  matter  of  fact,  the  approximation  (8.1)  corresponds  to  the  zero- 
Froude-number  limit  of  the  first-order  slender-ship  approximation  obtained  in 
Noblesse  (1978)  for  the  problem  of  the  wave  resistance  of  a  ship. 

The  more  complex  approximation  (8.2)  essentially  corresponds  to  a  second- 
order  approximation.  In  fact,  in  the  particular  case  of  longitudinal  translation 
of  a  slender  body,  we  have  |/^[f(x)-f (C)]Gnda| <<| f (t) | ,  and  the  approximation 
<P(S)  is  asymptotically  equivalent  to  the  straightforward  second-order  slender- 
body  approximation 

*(2)(t)  *  f(t)  +  [  [f (x)-f (C) ]G  (x,t)da(x)  (8.3) 

-'b  n 

given  by  the  recurrence  relation  (3.6),  with  (and  $^“f).  However,  the 

approximation  (8.2)  has  broader  applicability,  as  is  suggested  by  the  results  of 
calculations  for  longitudinal  and  transversal  translation  of  ogives  reported  in 
figures  9,  10,  and  11,  and  by  the  fact  that  the  potential  #(<t)  is  exact  in  the 
particular  case  of  translation  of  ellipsoids,  as  was  shown  previously.  Indeed, 
these  results  indicate  that  the  explicit  approximation  fl(t)  may  be  of  practical 
usefulness  for  a  large  class  of  bodies  and  body  motions,  excluding  however 
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transversal  motions  of  flat  disk-like  bodies.  More  precisely,  if  1,8,  and  6  are 
the  three  main  dimensions  of  a  body,  with  1>8>6 ,  it  was  predicted  that  the  approxi¬ 
mation  <P(t)  could  lead  to  numerical  difficulties  for  motion  in  the  direction  of 
the  minor  dimension  5  if  <5£( . 65- .  3f? ) 8 .  This  rough  tentative  criterion  for 

applicability  of  the  approximation  $(£)  does  not  exclude  usual  ship  hull  forms, 
for  instance  (typical  values  of  8  and  6  for  double-hull  forms  are  8“.15  and  6-.1). 

Simple  explicit  approximations  analogous  to  the  above  approximations  (8.1), 
(8.2),  and  (8.3)  can  likewise  be  obtained  for  the  analogous,  although  more  complex 
problems  of  potential  flow  about  bodies  in  the  presence  of  a  free  surface.  In 
particular,  first-and  second-order  slender-ship  approximations  corresponding  to 
approximations  (8.1)  and  (8.2),  respectively,  are  given  in  Noblesse  (1978)  for 
the  problem  of  the  wave  resistance  of  a  ship.  The  integral  equation  (3.2)  and 
the  explicit  approximation  (8.2)  have  also  been  extended  to  the  problem  of  poten¬ 
tial  flow  about  a  body  in  regular  water  waves  in  Noblesse  (1980). 
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Appendix  Integral  identities  for  the  interior  problem 

and  for  the  combined  exterior-interior  problems 

Although  the  present  study  is  primarily  concerned  with  the  "exterior 
potential-flow  problem",  that  is  the  problem  of  potential  flow  about  a 
body,  it  may  be  interesting  to  list  here  the  integral  identities  corresponding 
to  equations  (2.5a,b,c)  and  equation  (2.7)  for  the  "interior  potential",  i 
say,  defined  in  the  interior  domain  d^.  The  integral  identities  correspond¬ 
ing  to  equations  (2.5a,b,c)  take  the  forms 


=  f  GV 2<j>  1dv  +  (G^-^G  )d« 

*  Ja  Jk  n  n 


(Ala) 


for  £  in  d^-b,  that  is  inside  the  surface  b. 


0  *  I  GV  <J>  dv  +  (G$  -<J>  G  )da 

7  A  K  n  n 


(Alb) 


for  £  in  d-b,  i.e.  outside  b,  and 


li  2i  f  -ll 

^7  =  GV  <j>  dv  +  (Qp  -<P  G  )da 

2  *  Ja  Jk  n  n 


(Ale) 


for  £  exactly  on  (smooth)b.  The  integral  identity  corresponding  to  equation 
(2.7)  takes  the  form 


0  =  f  GV^dv  +  I  [Gd»±-(<|»1-<>»*)G  ]da  . 

J  A  Jk  n  n 


This  Integral  identity,  like  equation  (2.7),  is  valid  for  any  point  £,  either 
inside,  outside,  or  exactly  on  the  surface  b,  and  indeed  is  equivalent  to 
the  set  of  the  three  classical  identities  (Ala,b,c). 

If  we  add  the  integral  identities  (2.5a)  and  (Alb),  we  may  obtain  the 
relation 


=  GV2*dv  +  GV24> ^dv  +  |  [G (4» i-ij)  )+(<(>-<(ii)G  ]da. 

*  h  h  Jk  n  n  n 


Addition  of  the  integral  identities  (2.5b)  and  (Ala)  yields  the  same  relation, 
except  for  the  fact  that  <f>^  on  the  left  side  is  replaced  by  <j>^.  We  there¬ 
fore  have  the  relation 


=  GV2<f>dv  +  I  {G(4>±  -  <|>  )-J-(4>— 4> ±) G  ]da, 

Jh  n  n  n 


where  <j>  on  the  left  side  and  in  the  first  integral  on  the  right  side  clearly 
corresponds  to  <{>  or  4*1  for  points  outside  or  inside  the  surface  b,  respec¬ 
tively.  Naturally,  the  integral  relation  (A. 3)  can  also  be  obtained  by  add¬ 
ing  identities  (2.7)  and  (A2) .  Indeed,  this  yields 


(1-CjM*  +  C 


*  L 


GV  4>dv  + 


f  i® 

Jh  n  n 


)+(*-$  ^  ]da 


where  C  =-/,  G  da.  It  can  easily  be  seen  from  equations  (2.6)  and  (2.3)  that 
1  b  n  i 

the  expression  (l-C^H^+C^iJ)^  is  identical  to  or  <J>^  for  t,  outside  or  in¬ 
side  the  surface  b,  respectively,  so  that  the  above  relation  is  identical 
to  relation  (A3) . 

The  integral  relation  (A3)  expresses  the  potential  $0i)  1°  the  entire 

2 

space  d+d^  in  terms  of  a  volume  distribution  of  sources  with  density  V  $, 
and  surface  distributions  of  sources  and  normal  dipoles  on  the  surface  b, 
with  densities  <j>^~4>n  and  d>— <t> ^ >  respectively.  Two  classical  results  in  po¬ 
tential  theory  immediately  follow  from  relation  (A3),  namely:  (i)  a  distri¬ 
bution  of  normal  dipoles,  with  strength  6  say,  on  a  surface,  S  say,  generates 

a  potential  <t>  3  /_ 6G  da  whose  value  is  discontinuous  across  the  surface  S 
®  n  e  i 

(specifically,  we  have  $  -4>  mS ,  where  superscripts  e  and  i  refer  to  the  "exterior" 
and  "interior"  sides  of  S,  respectively;  the  "interior"  side  being  that  into 
which  the  unit  normal  vector  to  S  is  drawn),  and  (ii)  a  distribution  of  sources, 
with  strength  a  say,  over  a  surface  S  generates  a  potential  $  3  /gGoda 

whose  normal  derivative  $  is  discontinuous  across  S  (specifically,  we  have 

Ai  ,e  .  n 

$  —<p  *<r). 
n  n 


DISTRIBUTION  LIST  FOR  UNCLASSIFIED 
TECHNICAL  REPORTS  AND  REPRINTS  ISSUED  UNDER 
CONTRACT  N0Q014-78-C-0169  TASK  NR  062-525 


Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314  (12  copies) 

Professor  Bruce  Johnson 
U.S.  Naval  Academy 
Engineering  Department 
Annapolis,  MD  21402 

Library 

U.S.  Naval  Academy 
Annapolis,  MD  21402 

Technical  Library 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Annapolis  Laboratory 
Annapolis,  MD  21402 

Professor  C.  S.  Yih 
The  University  of  Michigan 
Department  of  Engineering  Mechanics 
Ann  Arbor,  MI  48109 

Professor  T.  Francis  Ogilvie 
The  University  of  Michigan 
Department  of  Naval  Architecture 
and  Marine  Engineering 
Ann  Arbor,  MI  48109 

Office  of  Naval  Research 
Code  211 

800  N.  Ouincy  Street 
Arlington,  VA  22217 

Office  of  Naval  Research 
Code  438 

800  N.  Quincy  Street 

Arlington,  VA  22217  (3  copies) 

Office  of  Naval  Research 
Code  473 

800  N.  Quincy  Street 
Arlington,  VA  22217 


NASA  Scientific  and  Technical 
Information  Facility 
P,  0.  Box  8757 

Baltimore /Washington  International 
Airport 

Maryland  21240 

Professor  Paul  M.  Naghdi 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  CA  94720 

Librarian 

University  of  California 
Department  of  Naval  Architecture 
Berkeley,  CA  94720 

Professor  John  V.  Wehausen 
University  of  California 
Department  of  Naval  Architecture 
Berkeley,  CA  94720 

Library 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Code  522.1 
Bethesda,  MD  20084 

Mr.  Justin  H.  McCarthy,  Jr. 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Code  1552 

Bethesda,  MD  20084 

Dr.  William  B.  Morgan 
David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Code  1540 

Bethesda,  MD  20084 

Director.  Office  of  Naval  Research  — 

Eas t em/Cent ral  Regional  Office  (Boston) 
Building  114,  Section  D 
666  Summer  Street 
Boston,  MA  02210 


-1- 


Library 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Technical  Library 
Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22418 

Technical  Documents  Center 

Army  Mobility  Equipment  Research  Center 

Building  315 

Fort  Belvoir,  VA  22060 

Technical  Library 

Webb  Institute  of  Naval  Architecture 
Glen  Cove,  NY  11542 

Dr.  J.  P.  Breslin 
Stevens  Institute  of  Technology 
Davidson  Laboratory 
Castle  Point  Station 
Hoboken,  NJ  07030 

Professor  Louis  Landweber 
The  University  of  Iowa 
Institute  of  Hydraulic  Research 
Iowa  City,  IA  52242 

Fenton  Kennedy  Document  Library 
The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  MD  20810 

Lorenz  G.  Straub  Library 

University  of  Minnesota 

St.  Anthony  Falls  Hydraulic  Laboratory 

Minneapolis,  MN  55414 

Library 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Technical  Library 

Naval  Underwater  Systems  Center 

Newport ,  RI  02840 

Engineering  Societies  Library 
345  East  47th  Street 
New  York,  NY  10017 


The  Society  of  Naval  Architects  and 
Marine  Engineers 

One  World  Trade  Center,  Suite  1369 
New  York,  NY  10048 

Technical  Library 

Naval  Coastal  Systems  Laboratory 

Panama  City,  FL  32401 

Professor  Theodore  Y.  Wu 
California  Institute  of  Technology 
Engineering  Science  Department 
Pasadena,  CA  91125 

Director,  Office  of  Naval  Research  - 
Western  Regional  Office  (Pasadena) 

1030  E.  Green  Street 
Pasadena,  CA  91101 

Technical  Library 
Naval  Ship  Engineering  Center 
Philadelphia  Division 
Philadelphia,  PA  19112 

Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709 
Editor 

Applied  Mechanics  Review 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  TX  78206 

Technical  Library 

Naval  Ocean  Systems  Center 

San  Diego,  CA  92152 

0NR  Scientific  Liaison  Group 
American  Embassy  -  Room  A-407 
AP0  San  Francisco,  CA  96503 

Librarian 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Defense  Research  and  Development  Attache 
Australian  Embassy 
1601  Massachusetts  Avenue,  NW 
Washington,  DC  20036 


Librarian  Station  5-2 
Coast  Guard  Headquarters 
NASSIF  Building 
400  Seventh  Street,  SW 
Washington,  DC  20591 

Library  of  Congress 

Science  and  Technology  Division 

Washington,  DC  20540 

Dr.  A.  L.  Slafkosky 
Scientific  Advisor 
Commandant  of  the  Marine  Corps 
Code  AX 

Washington,  DC  20380 

Maritime  Administration 
Office  of  Maritime  Technology 
14th  &  E  Streets,  NW 
Washington,  DC  20230 

Maritime  Administration 
Division  of  Naval  Architecture 
14th  &  E  Streets,  NW 
Washington,  DC  20230 

Dr.  G.  Kulin 

National  Bureau  of  Standards 
Mechanics  Section 
Washington,  DC  20234 


Professor  Justin  E.  Kerwin 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Cambridge,  MA  02139 

Professor  Phillip  Mandel 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Cambridge,  MA  02139 

Professor  J.  Nicholas  Newman 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Room  5-324A 
Cambridge,  MA  02139 

Professor  Francis  Noblesse 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Cambridge,  MA  02139 

Professor  Ronald  W.  Yeung 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Cambridge,  MA  02139 

Dr.  Robert  K.  C.  Chan 
JAYC0R 

1401  Camino  Del  Mar 
Del  Mar,  CA  92014 


Naval  Research  Laboratory 
Code  2627 

Washington,  DC  20375  (6  copies) 

Library 

Naval  Sea  Systems  Command 
Code  09GS 

Washington,  DC  20362 


Mr.  Marshall  P.  Tulin 
Hydronautics ,  Incorporated 
7210  Pindell  School  Road 
Laurel,  MD  20810 

Naval  Ship  Engineering  Center 
Code  6110 

Washington,  DC  20362 


Mr.  Thomas  E.  Peirce 
Naval  Sea  Systems  Command 
Code  03512 

Washington,  DC  20362 

Professor  Paul  Lleber 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  CA  94720 

Professor  C.  C.  Mei 

Massachusetts  Institute  of  Technology 
Department  of  Civil  Engineering 
Cambridge,  MA  02139 


Naval  Ship  Engineering  Center 
Code  6114 

Washington,  DC  20362 

Naval  Ship  Engineering  Center 
Code  6136 

Washington,  DC  20362 

Commandant,  U.S.  Coast  Guard 
G-0MI/31 


2100  Second  Street,  Southwest 
Washington,  D.C.  20593 


-3- 


